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1 Introduction 



High energy nuclear collisions will produce many thousands charged particles 
in the central rapidity interval usually covered by tracking detectors. Among 
the measured charged hadrons are those which are: 

• directly produced in the process of strong interaction between two nu- 
clei, 

• products of decay of resonances during the hadronization process and 

• products of weak decays of heavy hadrons. 

The characteristic distance between emission points of particles originat- 
ing from the first two processes is of the order 1 to 100 fm; this distance can 
be measured only in the statistical way using particle interferometry. How- 
ever, the characteristic distance between weak decay point and the collision 
region is often well measurable using recent experimental technics. 

In this note we discuss a possible method of statistical selection of weakly 
decaying particles in the experimental case when the characteristic decay 
distance is of the same order as experimental resolution of its measurement. 
In this case the main problem is to make efficient rejection of background 
due to the particles originating directly from the interaction. 

2 General Description 

We consider here two body decay of a given particle and assume that mo- 
menta and charges of both decay products are measured by the experimental 
set up. Typical examples are neutral strange or charm particle decays into 
charged particles (V topology) which can be easily measured. 
The basic points of the method are the following: 

• for each combination of two tracks coming from the same interaction, 
which are potential candidates for products of the decay the hypothet- 
ical decay point is fitted and the kinematical fit is done. 

• values of cut-variables are calculated for each pair and the probability 
density function (pdf) for accepted pairs from the same interaction in 
the cut-space (see below) is obtained ('raw signal' pdf). 
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• the 'raw background' pdf in the cut-space is calculated using combi- 
nations of tracks from different events, 

• the acceptance region in the cut-space is selected using 'raw back- 
ground' pdf, expected 'signal' and 'background' multiplicities and num- 
ber of collisions, 

• the 'signal' in the acceptance region is calculated as a difference between 
'raw signal' and 'raw background' pdfs multiplied by a 'raw signal' 
average multiplicity. 

The important element of the method is an introduction of the cut-space, 
using which transparent and well defined method for the selection of the 
best acceptance region (cuts) for the given experimental case is developed. 
The procedure should be used separately for different y — pr regions. The 
weighting factors, due to applied cuts in the cut-space necessary to obtain 
final results are trivial. 

3 Cut— variables 

In this Section we describe a set of five cut-variables which are later used 
to select the best acceptance region for separation of 'signal' from the 'back- 
ground'. The cut-variables are defined in a specific way, which grantees 
several important features of these variables calculated for the 'signal' (real 
V decays): 

• the variables are independent from each other, 

• the variables are independent of kinematical properties of the decaying 
V particle (eg. y or pt), 

• the 'signal' pdf in cut-variables is independent of an experimental res- 
olution, 

• the 'signal' pdf is uniform; the variables range between and 1. 

• the cuts in the cut-variables do not affect acceptance in the physically 
important variables eg. rapidity or transverse momentum. 
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In the following a specific set of 5 cut-variables is described. 



1. The cumulative of the decay vertex fit. 

The cumulative variable of the vertex fit is defined as: 



(1) 







where /(m, v) is chi-square probability density function for v degrees 
of freedom. The x^ value for each decay is calculated by decay vertex 
fit procedure. 

2. The cumulative of the kincmatical fit. 
The cumulative x^ variable is defined as: 



where fiu^v) is chi-square probability density function for u degrees 
of freedom. The x^ value for each decay is calculated using standard 
method of Least Squares estimation with constrains (kinematical fit). 
In of measured momenta vectors of both decay products and 

polar, and azimuthal, </)y, angles of the decay point (defined in the 
spherical system with the origin in the interaction point) the number of 
degrees of freedom is i/ — ?>. In the limit when there is no information 
on 9v and 0y = 1. 

3. The polar angle of the decay product. 
This cut-variable is defined as: 



where 9* is an angle between momentum vector of positive (or negative) 
decay product and V momentum vector calculated in the V center of 
mass system. The values obtained from the kinematical fit are used for 
the calculation. 




(2) 



Ce* =0.5(cosr + 1), 



(3) 
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4. The azimuthal angle of the decay product. 
This cut- variable is defined as: 




27r 



(4) 



where (j) is an angle between a projection of the beam direction and 
the projection of the momentum vector of the positive (or negative) 
decay product onto plane perpendicular to the V direction. The values 
obtained from the kinematical fit are used for calculations. 

5. The cumulative life time. 

The cumulative life time cut-variable is defined as: 





where i? is a distance between measured V decay point and the interac- 
tion point, r is a V life time in the V center of mass system and Pp^{R) 
is a probability density function of decay at the distance R calculated 
for a V with momentum pv- In the case of an ideal measurement of 
the decay distance R one gets: 



where Rq = pvTo/my and A is normalization factor. If the decay 
distance measurement is not an ideal one we get: 



where p{R, R') is a probability density that decay at a distance R! is 
measured as a decay at the distance R. The values obtained from 
kinematical fit are used for the calculations. 

The choice of the given set of cut-variables is not unique, especially vari- 
ous definitions of angular variables are possible. The final selection depends 
on the experimental situation. The variables should be defined in order to 
maximize nonuniformity of the 'raw background' pdf. 





(5) 



P,,{R) = PUR) = Ae-— = Ae-^, 



(6) 




(7) 
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4 Cut-space and Acceptance Region 



We define a 5-D cut-space as a set of points: 



C — (Cy, Ck-, Ce*, Cfj,, Cr), 



(8) 




(9) 



For the 'signal' pdf we obtain: 
Ps{C) = const{C) = 1 



(10) 



where Ps(C) is a 'signal' pdf at point C. 

The 'background' pdf, Pfc(C), is obviously nonuniform and therefore the 
ratio of the 'signal' to the 'background' depends on the point C or on the 
subspace selected in the cut-space. 

The average number of 'signal' events (averaged over all collisions) in the 
subspace of the volume V"''^'^ around any point C is given by: 



where < > is an average number of all 'signal' events, eg. all charged 
decays for which both decay products were measured. 

The corresponding number of 'background' events is given by: 



where < n;, > is an average number of 'background' events and p'^'^'^ is an 
average 'background' density in the subspace V°^. 

However, the average number of 'background' events and the 'background' 
pdf can not be directly measured due to unknown contamination of 'signal' 
events. The 'background' pdf, Pfe(C), can be safely estimated by 'raw back- 
ground' pdf, prb{C), obtained from combinations of tracks from different 
events providing that < Ug >«< rib > and that V decay products have 
similar y — pt distribution to that obtained for tracks originating from the 
main vertex. 

The only multiplicity which is measured is 'raw signal' multiplicity i.e. 
the average multiplicity of all track pairs, which are accepted as candidates 



< n",'" >= V"" < >, 



(11) 




(12) 
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for products of V decay. < nj-s >=< Ub > + < Us >. If the total 'signal' 
multiplicity is much smaller than total 'background' multiplicity one can 
safely write: 

< rib >~< rirs > ■ (13) 

The above approximation allows to calculate 'signal' multiplicity in the given 
acceptance (V"'^'^): 

< n^/'' >=< <f > - < >^< <f > -p^b^"^" < rirs >, (14) 

providing that the uncertainty of the < n"'^'^ > due to the used approximation 
(Eq. 13) is much smaller than < nf^ > i.e.: 

< n > 

< »> < > (15) 

<nb> 

or consequntly 

« 1- (16) 

Therefore in the region of cut-space in which the last inequahty is satisfied 
the 'signal' multiplicity can be expressed as: 

<ns>=<nrs>{rrT-X'tl- (17) 

For a given number of collisions, Nev, the 'signal' can be extracted from 
the 'background' in the subspace V"''^'^ when the statistical uncertainty of the 
< > is much smaller than < >. Assuming that < rirs >, VrT Prb^ 
are independent one gets: 

a{< n, >) = [(a(< n,, >)(jfrT-rrtnf+i< n.. > a(p,T))'+(< rirs > aiXnTY^'m 

Obviously the second term dominates in the above expression (for cen- 
tral collisions < nrs > is approximately fixed and therefore cr(< rirs >) is 
relatively small even for small number of collisions; number of background 
combinations increarses like N^^ whereas number of signal like Ngy and there- 
fore cr(j)":^'^) is \/Nev times smaller than a{p^). Neglecting the first and the 
last terms in the Eq. 18 we get: 

a{< ris >) =< rirs > 'ji'KT)- (19) 
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Let us try to estimate a{p'^f'). From definition we have: 
^cc ^ ^20) 

Due to the fact that v^^^ « rirs we can assume Poissonian flactuations of 
n"^'^ and therefore we get: 



acc 



^ acc ^ (21) 

Substituting n"^"^ = n^-s 'p^^ one gets: 



-^acc Y 

Therefore the minimum number of events neccesary for signal-background 
separation is given by: 

TV,, » ^ \ (23) 
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